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Abstract: The power Bonferroni mean (PBM) operator is a hybrid structure and can take the 
advantage of a power average (PA) operator, which can reduce the impact of inappropriate data given 
by the prejudiced decision makers (DMs) and Bonferroni mean (BM) operator, which can take into 
account the correlation between two attributes. In recent years, many researchers have extended the 
PBM operator to handle fuzzy information. The Dombi operations of T-conorm (TCN) and T-norm 
(TN), proposed by Dombi, have the supremacy of outstanding flexibility with general parameters. 
However, in the existing literature, PBM and the Dombi operations have not been combined for 
the above advantages for interval-neutrosophic sets (INSs). In this article, we first define some 
operational laws for interval neutrosophic numbers (INNs) based on Dombi TN and TCN and 
discuss several desirable properties of these operational rules. Secondly, we extend the PBM operator 
based on Dombi operations to develop an interval-neutrosophic Dombi PBM (INDPBM) operator, an 
interval-neutrosophic weighted Dombi PBM (INWDPBM) operator, an interval-neutrosophic Dombi 
power geometric Bonferroni mean (INDPGBM) operator and an interval-neutrosophic weighted 
Dombi power geometric Bonferroni mean (INWDPGBM) operator, and discuss several properties of 
these aggregation operators. Then we develop a multi-attribute decision-making (MADM) method, 
based on these proposed aggregation operators, to deal with interval neutrosophic (IN) information. 
Lastly, an illustrative example is provided to show the usefulness and realism of the proposed MADM 
method. The developed aggregation operators are very practical for solving MADM problems, as it 
considers the interaction among two input arguments and removes the influence of awkward data 
in the decision-making process at the same time. The other advantage of the proposed aggregation 
operators is that they are flexible due to general parameter. 


Keywords: interval neutrosophic sets; Bonferroni mean; power operator; multi-attribute decision 
making (MADM) 


1. Introduction 


While dealing with any real world problems, a decision maker (DM) often feels discomfort when 
expressing his\her evaluation information by utilizing a single real number in multi-attribute decision 
making (MADM) or multi-attribute group decision making (MAGDM) problems due to the intellectual 
fuzziness of DMs. For this cause, Zadeh [1] developed fuzzy sets (FSs), which are assigned by a 


Symmetry 2018, 10, 459; doi:10.3390/sym10100459 www.mdpi.com/journal/symmetry 


Symmetry 2018, 10, 459 2 of 32 


truth-membership degree (TMD) in [0,1] and are a better tool to present fuzzy information for handling 
MADM or MAGDM problems. After the introduction of F5s, different fuzzy modelling approaches 
were developed to deal with uncertainty in various fields [2-4]. However, in some situations, it 
is difficult to express truth-membership degree with an exact number. In order to overcome this 
defect and to express TMD in a more appropriate way, Turksen [5] developed interval valued FSs 
(IVFSs), in which TMD is represented by interval numbers instead of exact numbers. Since only TMD 
was considered in FSs or IVFSs and the falsity-membership degree (FMD) came automatically by 
subtracting TMD from one, it is hard to explain some complicated fuzzy information, for example, for 
the selection of Dean of a faculty, if the results received from five professors are in favor, two are against 
and three are neither in favor nor against. Then, this type of information cannot be expressed by FSs. So, 
in order to handle such types of information, Atanassov [6] developed intuitionistic fuzzy sets (IFSs), 
which were assigned by TMD and FMD. Atanassov et al. [7] further enlarged IFSs and developed 
the interval valued IFS (IVIFSs). However, the shortcoming of FSs, IVFSs, IFSs and IVIFSs are that 
they cannot deal with unreliable or indefinite information. To solve such problems, Smarandache [8,9] 
developed neutrosophic sets (NSs). In neutrosophic set, every member of the domain set has TMD, 
an indeterminacy-membership degree (IMD) and FMD, which capture values in JO”, 1*[. Due to the 
containment of subsets of ]0~, 1*[ in NS, it is hard to utilize NS in real world and engineering problems. 
To make N&Ss helpful in these cases, some authors developed subclasses of NSs, such as single valued 
neutrosophic sets (SVNSs) [10], interval neutrosophic sets (INSs) [11,12], simplified neutrosophic sets 
(SNSs) [13,14] and so forth. In recent years, INSs have gained much attention from the researchers 
and a great number of achievement have been made, such as distance measures [15-17], entropies of 
INS [18-20], correlation coefficient [21-23]. The theory of NSs has been extensively utilized to handle 
MADM and MAGDM problems. 

For the last many years, information aggregation operators [24-27] have stimulated much 
awareness of authors and have become very dominant research topic of MADM and MAGDM 
problems. The conventional aggregation operators (AGOs) proposed by Xu, Xu and Yager [28,29] can 
only aggregate a group of real numbers into a single real number. Now these conventional AGOs were 
further extended by many authors, for example, Sun et al. [30] proposed the interval neutrosophic 
number Choquet integral operator for MADM and Liu et al. [31] developed prioritized ordered 
weighted AGOs for INSs and applied them to MADM. In addition, some decision-making methods 
were also developed for MADM problems, for example, Mukhametzyanov et al. [32] developed a 
statistically based model for sensitivity analysis in MADM problems. Petrovic et al. [33] developed 
a model for the selection of aircrafts based on decision making trial and evaluation laboratory and 
analytic hierarchy process (DEMATEL-AHP). Roy et al. [34] proposed a rough relational DEMATEL 
model to analyze the key success factor of hospital quality. Sarkar et al. [35] developed an optimization 
technique for national income determination model with stability analysis of differential equation in 
discrete and continuous process under uncertain environment. These methods can only give a ranking 
result, however, AGOs can not only give the ranking result, but also give the comprehensive value of 
each alternative by aggregating its attribute values. 

It is obvious that, different aggregation operators have distinct functions, a few of them can reduce 
the impact of some awkward data produced by predispose DMs, such as power average (PA) operator 
proposed by Yager [36]. The PA operator can aggregate the input data by designating the weight 
vector based on the support degree among the input arguments, and can attain this function. Now the 
PA operator was further extended by many researchers into different environments. Liu et al. [37] 
proposed some generalized PA operator for INNs, and applied them to MADM. Consequently, some 
aggregation operators can include the interrelationship between the aggregating parameters, such as 
the Bonferroni mean (BM) operators developed by Bonferroni [38], the Heronian mean (HM) operator 
introduced by Sykora [39], Muirhead Mean (MM) operator [40], Maclaurin symmetric mean [41] 
operators. In addition, these aggregation operators have also been extended by many authors to deal 
with fuzzy information [42-46]. 
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For aggregating INNs, some AGOs are developed by utilizing different T-norms (TNs) and 
T-conorms (TCNs), such as algebraic, Einstein and Hamacher. Usually, the Archimedean TN and TCN 
are the generalizations of various TNs and TCNs such as algebraic, Einstein, Hamacher, Frank, and 
Dombi [47] TNs and TCNs. Dombi TN and TCN have the characteristics of general TN and TCN by a 
general parameter, and this can make the aggregation process more flexible. Recently, several authors 
defined some operational laws for IFSs [48], SVNSs [49], hesitant fuzzy sets (HFSs) [50,51] based on 
Dombi TN and TCN. In practical decision making, we generally need to consider interrelationship 
among attributes and eliminated the influence of awkward data. For this purpose, some researchers 
combined BM and PA operators to propose some PBM operators and extended them to various 
fields [52-55]. The PBM operators have two characteristics. Firstly, it can consider the interaction 
among two input arguments by BM operator, and secondly, it can remove the effect of awkward data 
by PA operator. The Dombi TN and TCN have a general parameter, which makes the decision-making 
process more flexible. From the existing literatures, we know that PBM operators are combined with 
algebraic operations to aggregate IFNs, or IVIFNs, and there is no research on combining PBM operator 
with Dombi operations to aggregate INNs. 

In a word, by considering the following advantages. (1) Since INSs are the more précised class by 
which one can handle the vague information in a more accurate way when compared with FSs and all 
other extensions like IVFSs, IFSs, [VIFSs and so forth, they are more suitable to describe the attributes 
of MADM problems, so in this study, we will select the INSs as information expression; (2) Dombi TN 
and TCN are more flexible in the decision making process due to general parameter which is regarded 
as decision makers’ risk attitude; (3) The PBM operators have the properties of considering interaction 
between two input arguments and vanishes the effect of awkward data at the same time. Hence, the 
purpose and motivation are that we try to combine these three concepts to take the above defined 
advantages and proposed some new powerful tools to aggregate INNs. (1) we define some Dombi 
operational laws for INNs; (2) we propose some new PBM aggregation operators based on these new 
operational laws; (3) we develop a novel MADM based on these developed aggregation operators. 

The following sections of this article are shown as follows. In Section 2, we review some basic 
concepts of INSs, PA operators, BM operators, and GBM operators. In Section 3, we review basic 
concept of Dombi TN and TCN. After that, we propose some Dombi operations for INNs, and discuss 
some properties. In Section 4, we define INDPBM operator, INWDPBM operator, INDPGBM operator 
and INWDPGBM operator and discuss their properties. In Section 5, we propose a MADM method 
based on the proposed aggregation operators with INNs. In Section 6, we use an illustrative example 
to show the effectiveness of the proposed MADM method. The conclusion is discussed in Section 7. 


2. Preliminaries 


In this part, some basic definitions, properties about INSs, BM operators and PA operators 
are discussed. 


2.1. The INSs and Their Operational Laws 


Definition 1. Let O be the domain set [8,9], with a non-specific member in QO. expressed by 0. A NS NS in Qis 
expressed by 


NS = {(3,tx<(®),ixg(®), feg(®) fo € OF, (1) 


where, t=(0),ix=(0) and fo. (2) respectively express the TMD, IMD and FMD of the element 0 € U to 


the set NS. For each point 3 € U, we have, to (2), I< (°) fr (3) € ]O-,1+[and0- < t=<(0) bi=-(0) 4 


fyg(2) <3". 


The NS was predominantly developed from philosophical perspective, and it is hard to be applied 
to engineering problems due to the containment of subsets of |0~, 1* |. So, in order to use it more easily 
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in real life or engineering problem, Wang et al. [8] presented a subclass of NS by changing ]0~,1*[ to 
[0,1] and was named SVNS, and is defined as follow: 

Definition 2. Let O be the domain set [10], with a non-specific member in O. expressed by 0. A SVNS SV inQ 
is expressed by 


SV = { (ii, te), ig¢®), fen) ) [8 € O}, (2) 


where t=>(0),ig;(0) and faz) express the TMD, IMD and FMD of the element 0 € Oto the se t SV 
respectively. For each point 0 € O, we have, t=;(0),iz;(), faz(®), € [0,1] and 0 < t=5(2) tig (0) 4 
fg(2) <3 


In order to define more complex information, Wang et al. [9] further developed INS which is 
define as follows: 


Definition 3. Let Q be the domain set and € Q [11]. Then an INS IN in Qis expressed by 


IN = { (3, TRyy (8), ID py (8), Flin (2) joe Oh, (3) 


where, TR; (0), ID;y (0) and FL;,,(0) respectively, express the TMD, IMD and FMD of the element 
0 € Q to the set IN. For each point 0 € U, we have, TRin(®),IDin(0), FLin(®) © [0,1] and 
0 < maxID ny (0) + maxID ny (0) + maxFL yy (0) < 3. 


— L u L —u Lu 
For computational simplicity, we can use in = (FR ,ITR hb {1D ,ID | [Fr ,FL }) to 
express an element in inan INS, and the element in is called an interval neutrosophic number (INN). 


Where TR TR | C [0,1], {TD ID | C [0,1], [Fr /FL | C [0,ljandO<TR +ID +FL <3. 


L—u L —u L —u _ 
Definition 4. Let in, = ( TR, TR, [ Dy, ID, h [Foy FL; | ) and ing a 


L uU L u L U 
TR z,TR, |, |ID2,ID> |, Fz, Fl | ) be any two INNs [12], and € > 0. Then the operational 
laws of INNs can be defined as follows: 


: : L L LL u u u—u LL —u—u LL —Uu—u 
(1) iny ® ing = ( | TR, + TR» — TR, TRy, TR, +TR> — TR, TR |, | ID, ID, ID, ID, |, |FL,FLy,FL, FL, |); (4) 


- - L——L u—u L L L—=L —u u u. 
(2) in, @ ing = (FR: TR», TR; TR» | [i + ID) —1D,ID;,1D; +ID, — 1D, iD: |, [Fr if FL, AL. FL, 4 FE, 6 er, | (5) 


om) (rR) f(y a(t] = (mE) (8 1) m 
1 (menor) [my wey} [ey ey]). 


L. u L Uu L U 
Definition 5. Let in = ( TR ,ITR | {TD ,ID | Fr ,FL | ) [42], be an INN. Then the score function 


S (in) and accuracy function A (in) can be defined as follows: 


+1 Et ; (8) 


in) = +TR ID +ID FL +FL 
2 2 2 
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+ J 5 7 : (9) 


—L =u L u L u 
in) ==" +TR ID +I1D FL +FL 


In order to compare two INNs, the comparison rules were defined by Liu et al. [36], which can be 
stated as follows. 


= L =u L =u L =u = 

Definition 6. Let in, = ( TR, TR, | 'TDy, ID, | [Foy FL; | ) and ing a 
L 

( TR, TR |, 


| 
(1) eo, 
fe 
Ve 


L U 
iD. ID, >|; FL, FLy | ) be any two INNs [42]. Then we have: 


(inz), then in is better than ino, and denoted by i in, > ind; 
If Si (im 


S 
S (inz), and A (im) >A (in2), then iny is better than ino, and denoted by iny > ino; 
ifs (im S 


(inz), and A (im) =A (in2), then iny is equal to in, and denoted by iny = ino. 


: L —u L —u L —u = 
Definition 7. Let in, = ( TR, TR, i TD, ID, | [Foy FL; | ) and ing = 


L u L u L U 
(| FR TR, 1B, TB, |, [Fl:,FT2 | ) be any two INNs [15]. Then the normalized Hamming 


distance between ny and np is described as follows. 


D(im,, in) = z(|FR TR, 


U u 
| TR, SIR, 


L L 
iD TD> 


u u 
iD 1D, 


L L 
Fr, EL 


U U 
Fr, FL, 


) ao) 


2.2. The PA Operator 


The PA operator was first presented by Yager [36] and it is described as follows. 


Definition 8. For positive real numbers ¢,(h = 1,2,...,1) [36], the PA operator is described as 


l 
Py (1+ T(@n))@n 
PA( 91, 92,-+ +, 1) = = , (11) 


l 
where, T(9,) = YL sup(@n, Gy), and sup(n, Py) is the degree to which yy, supports gy. The support 


yi 


degree (SPD) satisfies the following properties. 


(1) sup(~n, Gy) = sup(—y, Pn); 
(2) sup(@p, @y) € [0,1]; 
(3) sup(n, (y) > sup(%c, a), if | on — %y| < | fc — goal. 


2.3. The BM Operator 


The BM operator was initially presented by Bonferroni [38], and it was explained as follows: 


Definition 9. For non-negative real numbers gy,(h = 1,2,...,1), and x,y > 0 [38], the BM operator is 
described as 


1 1 x+y 
yy cist) | (12) 


BM*Y (01, 92,--+, 91) = ( =] 
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The BM operator ignores the importance degree of each input argument, which can be given by 
decision makers according to their interest. To overcome this shortcoming of BM operator, He et al. [52] 
defined the weighted Bonferroni mean (WBM) operators which can be explained as follows: 


Definition 10. For positive real numbers g,(h = 1,2,...,1) and x,y > 0 [52], then the weighted BM operator 
(WBM) is described as 


1 


ee x+y 
Ky Ks xy 1 
1 = Kp Ph, #) , ( 3) 


l l 


1 
WBM™ (01, 92,.-+, 1) = (ey » 
h=1s=1,h4s 


where % = (K1,K,...,%,)' is the importance degree of every o,(h = 1,2,...,1). 


The WBM operator has the following characteristics: 


T 
Theorem 1. (Reducibility) If the weight vector is K = (}. i, ae t) , then 


h=1s=1,z4s 


1 1 x+y 
WBM*" (91, 92,---,1) = (2, ae vist) (14) 
= BM*Y (91, 92,.-+,m). 


Theorem 2. (Idempotency) Let oy, = 9, (h = 1,2,...,1). Then BM*4 (91, 92,..., 91) = 9. 
Theorem 3. (Permutation) Let ( (1, 92,.-., 01) be any permutation of (Z1',Z2',...,Z)'). Then 


WMB*4 (Zy',Zp',...,Z)') = WBM(1, 92,---, 91): (1) 


Theorem 4. (Monotonicity) Let o;, > Kj,/(h =1,2,...,1). Then 


WBM*” (01, 02,---, 1) > WBM* (Ky’, Ko! yeaa Ky’). (16) 


Theorem 5. (Boundedness) The WBM*¥ lies in the min and max operators, that is, 


min((1, PI2r- +67 (21) < WBM*” (04, fI2r- +07 21) < max(91, §I2r-++7 (21): (17) 


Similar to BM operator, the geometric BM operator also considers the correlation among the input 
arguments. It can be explained as follows: 


Definition 11. For positive real numbers o,(h = 1,2,...,1) and x,y > 0 [53], the geometric BM operator 
(GBM) is described as 


i i 1 
IL I] @entyps)?. (18) 


GBM*" (1, (02,---, 91) = : Ty 
h=1s=1,h4s 


The GBM operator ignores the importance degree of each input argument, which can be given 
by decision makers according to their interest. In a similar way to WBM, the weighted geometric 
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BM (WGBM) operator was also presented. The extension process is same as that of WBM, so it is 
omitted here. 


The definition of power Bonferroni mean (PBM) and power geometric Bonferroni mean (PGBM) 
operators are given in Appendix A. 


3. Some Operations of INSs Based on Dombi TN and TCN 
Dombi TN and TCN 


Dombi operations consist of the Dombi sum and Dombi product. 


Definition 12. Let S and & be any two real numbers [47]. Then the Dombi TN and TCN among S and ¥ are 
explained as follows: 


Ta(@,8) = re (19) 


Tp*(S,X) =1 ~, (20) 
T 


where, 1 > 1, and (3,%) € [0,1] x [0,1]. 


According to the Dombi TN and TCN, we develop a few operational rules for INNs. 


L U L U L uU = 
Definition 13. Let in = (TR ,ITR |, | ,ID |,|FL ,FL ) in = 


L uU L uU 
( [FR TR; ,|ID,,1D, [ 


L =u = L =u L =u L —u 
FL,,FL, | ) and ing = ( [FR TR: |, 1B, 1B, FL, FL; ) 
be any three INNs and ® > 0. Then, based on Dombi TN and TCN, the following operational laws are developed 
for INNs. 


= 1 1 1 
int ® in2 = i= A | , 


1, NP ==L \ \y =u \ =u )\)\y 
TR2 TR TR2 
(1) I+ ZL ot: =r Tet: = gs 7 
1-TR, 1-TR2 1-TR, 1-TR2 


iny ® ing = ( I T I 
asl \T ETN —u\T =sU\7\ 7 aT, VY =t. \I\ F 
4 1-TRy p( 12 TRe 4 1-TRy , { 1TRo 14 1D} 4( 722 
(2) TR, TRE TRI TR 11D) 1-105 (22) 
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a\h = yr avi Ne aa 
(3) (Re) ) (Ce) LLG) ) CGE) ) Lay) 
ene ee } 
(CHE) ) “CCH)) 
= ( 1 a 1 . 1 1 wi 1 : 
() (CE) ) MCG) TL Gh) 0G) )) ay 


Now, based on these new operational laws for INNs, we develop some aggregation operators to 
aggregate IN information in the preceding sections. 


4. The INPBM Operator Based on Dombi TN and Dombi TCN 


In this part, based on the Dombi operational laws for INNs, we combine PA operator and BM 
to introduce interval neutrosophic Dombi power Bonferroni mean (INDPBM), interval neutrosophic 
weighted Dombi power Bonferroni mean, interval neutrosophic Dombi power geometric Bonferroni 
mean (INDPGBM) and interval neutrosophic weighted Dombi power Bonferroni mean (INWDPGBM) 
operators and discuss some related properties. 


4.1. The INDPBM Operator and INWDPBM Operator 


L =u L—u — 
Definition 14. Let in; = (FR , TR; | 5; ,ID; | Fa; , FL; » (i =1,2,...,1), be a group of INNs, 
and x,y > 0. If 


x 
. 


= _ y xy 
AG) =| os Kei) i 05) 
1 ® (1+T(in,)) © (1+T(im,)) 
u=1 u=1 


then INDPBM*¥ is said to be IN Dombi power Bonferroni mean (INDPBM) operator, where T(in:) = 


wilt = = 1 : 
INDPBM* (in, ing, ...,im)) ~~) RP -] : 


if] 


1 — — — — — — 
@ Sup (inz, ins). Sup (inz, ins) is the support degree for inz from ings, which satisfies the following 
s=1sAz 

axioms: (1) Sup (inz,ins) € [0,1]; (2) Sup (inz,ins) = Sup(ins,inz); (3) Sup (inz,ins) > Sup(itta,iny)), 
if D inz, ins <D (ita, inp) ,inwhich D (ita, inp) is the distance measure between INNs ing and inp defined 
in Definition 7. 


In order to simplify Equation (25), we can give 


Ke) i 


2(1+7(m)) 
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l 
and call A = (Aj, Ao,... ,Ai)* is the power weight vector (PWV), such that A, > 0, @ A, = 1. It 
z=1 


turns Equation (25) into the following form 


INDPBM* (inj, img,...,im) = — & (JAiin;) @p (1Ajin;)" (7) 


= L =u L =u L —==L 
Theorem 6. Let in; = (FR , TR; |, 1B; TB; | FE; FL; | ),( = 1,2,...,1) be a group of INNs. 
Then the value obtained by utilizing Equation (25) is expressed as 


/ 1 | \ 
a( | 4 [ibe y yf ee d y fisfote yf ly a *__¥4 » 
x+y / =] = = VY +y 7A = u 
yf VIA we in| 2B YE | iy) 78 tn, | 2B 
1-TR 1-7R, / / 1-TR 1-TR; 
/ 1 L 
| / / 
/ / | 
/ / i 
i= af |iafFE fi = y io fiaf{tete af ds i/ as a 
7 a zy / » / =L\" - ==L\" an : x+ te / p> / =u \" ‘ =u \! i 
fled lan a om eB J ie / lin ea | + a 
/ ih i = i = | if i =U ‘ =U 
/ / ID; ID; / / ID; ID; (28) 
1 I 
P-l ae ay x y Pop ed x y 
I- Y \ie}—x y FY y —+—=|]] |I- If fl+]—x I] \ > 
tty of jaa | al aa xty of daa | ay aay 
sf Lan 1-FLi 1-FL Sf a | 1-FLi 1-FL 
| | IA; — jf, IA, ——f; | i | IA, = IA, ty 
| FLi FL; | | FLi FLj 


Proof. Proof of Theorem 6 is given in Appendix B. 


In order to determine the PWV A, we firstly need to determine the support degree among INNs. 
In general, the similarity measure among INNs can replace the support degree among INNs. That is, 


Sup (ini, im) =1- D(ini,inm) (i,m =1,2,...,1). (29) 


Example 1. Let in; = ((0.3,0.7], (0.2,0.4], [0.3,0.5]),inz = ({0.4,0.6], [0.1,0.3], [0.2,0.4]) and ing = 
([0.1,0.3], [0.4,0.6], [0.2,0.4]) be any three INNs, x = 1y = 1,y = 3. Then by Theorem 6 
in Equation (28), we can aggregate these three INNs and generate the comprehensive value in = 


L U L U L uU 
(FR ,ITR | i ,ID | Fz ,FL |) which is calculated as follows: 
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Step1. Determine the supports Sup(in:, inj) i,j = 1,2,3 by using Equation (29), and then 
we get Sup (in, ing) = Sup (ing, im; ) = 0.9, Sup (im, ins ) = Sup (ins, in) = 
0.933, Sup (ina, ins ) = Sup (ins, in ) = 
—- 3 —_—- —— 
Step2. Determine the PWV Because T(in:) = y Sup(inz, ins), and 
s=1sAz 


Sup(in ins] = Sup ins.in = 0.9, Sup im.ins} = Sup insin) = 0.933, up| inns} = Sup ins.ins} =1, 


_ Ay = = (7(m) ag !) — = 0.3269 
(r() +1) + (r(He) +1) + (Gm) 2) 
Ag = — (7(im2) = 1) — = 0.3346 
(F(R) 1) + (re) 41) +H) ay 
A3 = (7(7s) = !) = = 0.3385. 

(FC) 2) + (rH) -4) + RH) 9) 


Step 3. Determine the comprehensive value in = 


a 
( re | 


L U E u L u 
TR ,TR I {1D ,ID i [Fz EL |) by using 
Equation (28), we have 


3 , ; 
x) vm z+ a el= 4 } =+ 2 z|+ | : z+ a =|+ J A z+ 2 : 
lh TR TR, TR, TR: TR TR TR: TR 
% 3A,] | 3A,] 3A, a} 3A.) 3A, | 3As] —p 3A,] =| 34) 
1-TR; 1-TR; 1-TR, 1-TR2 1-TR 1-TR; 1-TR2 1-TR 
+ | : zt 2 z|+ | u at 2 z|+ | ! z+ 2 = |= 0.1281 
=r _ 1 = = 1 
3A, TR 3A, se 3A, ERS. 3A, TRE 3A, HS. 3A, ee 
1-TR2 1-TR3 1-TR; 1-TRi 1-TR3 1-TR2 


[oo 


32-3 : 1 1 
1/14 cr 1/ yy 1 ——; 4 —y = 0.2590 
ipa 3A ( a 24,( 7 | 
iXj 1—-TR; 1-TR| 


Similarly, we can get n = ([0.2590, 0.5525], (0.2221, 0.4373], [0.2334, 0.4365]). 


= L =u L =u L ==L 
Theorem 7. (Idempotency) Let inj = (FR , LR; | 5; ,ID; | FE, ,FL; i Gi= 


1,2,...,1), be a group of INNs, if all inj(i =1,2,...,1) are equal, that is ins = in = 
18 Uu L u L uU 
[TR ,TR |B ,1D |. [FL ,FL ))@=12,...D, Then 


INDPBM*Y (im, ily, .. . in) = in. (30) 
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Proof. Proof of Theorem 7 is given in Appendix C. 


Theorem 8. (Commutativity) Assume that in’ is any permutation of iny(u = 1,2,...,1), then 


INDPBM* (in, in'y,... in!) — INDPBM* (iny,ino, os /in,) 


Proof. From Definition 14, we have 


a 
INDPBM*# (in, in’, ie) = — ¥ (1A'iin’;)” Qp (1A/jin'))” ; 
pial 
coy 
and 1 
xy 
INDPBM*# (iny,ins, 7 .,in}) = 7 3 (JAiin;)" Qp (JAjin;)" 
,j=1 
Tei 
Because, 
ie (1A'jin',)" SD (1A' jin’) = x (iain) Sp (1Ajinj)", 
,j=l rey 
T4j une 


Hence, INDPBM* (ins, in’, 7 . in’) = INDPBM*Y (ity, ing,. sia) 


— L u L u LL 
Theorem 9. (Boundedness) Let in; = (FRR i [1D;, 1D; i [FE,FE;| ).( ee ee) ee 


be a group of INNs, and in = max | TR; , TR; min ID; , ID; min FL; , FL; ,in 
1= = 


i=1 
T [a= =u 1 [==l =u 1 f==_L =u 
(min TR , TR; | mt ID; , 1D; | mi Fr; , FL; be Then 
1= {= — 


in < INDPBM (in, iny,...,in;) <in . 


Proof. Proof of Theorem 9 is given in Appendix D. 


(31) 


(32) 


Now, we shall study a few special cases of the INDPBM*Y with respect to x and y. (1) When 


y 0,7 > 0, then we can get 
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aa a 

/ 7 7 
P-l t x P-l t x 
1- 1+ 1 1 ———_—_—_——- 3 1 1 1 1 ——_—__—- ‘ 
(SV | =a (Se yf |S Sy (33) 
ij 1A, cies i ij 1A, = 
ID; | ID; 
/ 


(2) When x = 1,y > 0,7 > 0, then we can get 


~ = 
Y 
=li 1 1+(P-1 1 1 — 1 14+(P-1 1 1 — 
=e piel Ta 2 = SPU = ; 


1 1) 
Y y 
(34) 
F L 1 e : 1 
1- Y |i+(P-ix Yo |X y¥ }———[]]] }] ¥ [(P-e ¥ |Y yy |—.\ | , 
ijl _ Tn. \ ij-l _TD 
DP |B # || 
Di 1D; 
1 1 
Y x 
f 1 1 : 1 
i ay (Poke Yay | tll a eae |e 
i _FL i _FL 
oD |p ( FE oP | nF 
FLi FLi 


(3) When x = y = 1,y > 0, then we can get 


He, 
1 


INDPBM“° (im,in sits in) 
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/ im | Te) a aie / al f mm \ = i 
| - / IA, a IA; _ | me / IA, = IA; = 
/ / 1-TRi 1-TR; | / 1-TRi 1-TR; 


P-1 Jie] * , PHI / be J 7 
I- Y fisf—ex y 2 y + 7——||] fle yo fat ee yy + _, P (35) 


In the INDPBM operator, we can only take the correlation among the input arguments and cannot 
consider the importance degree of input arguments. In what follows, the INWPDBM operator shall be 
proposed to overcome the shortcoming of the INDPBM operator. 


L =u L =u L ==L 
Definition 15. Let in; = ( TR;,TR; I 75; ,1D; I Fr; , FL; | ) (i = 1,2,...,1), bea group of INNs, 
then the INWDPBM operator is defined as 


x 


INWDPBM™ (in ,ins....in} = P pa — in. | ®, 
SO ie (rin. +} 
zal 


bv, (7G) +1) - 


>. (rin. +} 7 


= l —- — 
where, T (ini) = Yo Sup (ini, inj), X,Y > 0,w = (w1,W2,... ,w)" is the importance degree of the INNs, 
J=liFj 


l 
such thatO <w, <1(z=1,2,...,l) and Y wp = 1. 
k=] 


= L —u L —u 
Theorem 10. Let inj = ( TR; , TR; | 1B; ,1D; | 


L ==L 
FL; , FL; | : (i = 1,2,...,1) be a group of INNs. 
Then the value obtained using Definition 15, is represented by 


INDPBM™" (im, ins,....im) 
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| 
| / 
f i ( yy : | (St i ya é + d 
\ oe Bf | (rein) La Jz (rn) +1) Ta mY ’ rey fli | tw,(Gn) +1) mY (temas) ( at 
\ / / — / oe — — _ 
\ / [ Yue lr (rv in. in. +0) 1-TR, (rin. +1 ue TR; } / [ / Yw.(rie. +) | Ye v.(Tin +1 iu- Tr, | 
\ / / a J) / { Vel 
| all 
1- ; (S i/ y Vf: # ¥ : / (S / y | ae + =! (37) 
aty a by, (rin) +) = _pyY 3 (tiny +1 it #|| sty fle Pty i( rn) +1) [8 py) ™ (Pin) y+) ‘Z D | 
/ rin. +))\\ lz (Pin 41 me (rin, in. +)\_ 1D wo } Ye. (rn. in. +)\_ 1D D, } 
/ % 
| nwa 
/ t 4 yy , 
/ / / { / 
(e- / 2 x y (E ee ae ee | x y 
i Yo fife 4 \ = Joy fufotx y V = = |. 
/ ere z | Iw (Gn) +1) (zy h (rcn)+1) [ze "ley / 2 / tw (7G) +1) [z le : w, [Pin + ) (+ i ] | 
| Ye. (rd. +0) FL Era FL, / / / Yw. (rd. +0) FL Ym (ri is in. +0) a | 
/ al = ) / / | = || 
N J 


Proof. Proof of Theorem 10 is similar to Theorem 6. 


Similar to the INDPBM operator, the INWDPBM operator has the properties of boundedness, 
idempotency and commutativity. 
4.2. The INDPGBM Operator and INWDPGBM Operator 

In this subpart, we develop INDPGBM and INWDPGBM operators. 


L =u L =U L ==L 
Definition 16. Let inj = ( |TR;,TR; | i ,1D; | Fr; FL; | ) (i = 1,2,...,1) be a group of INNs. 
Then the INDPGBM operator is defined as 
te 
- - 2-1 
\(TGnj)+1) U(T (fa) +1) 
a = it l =} (nie)+4) = (T(inz)+1) 
INDPGBM” (ity, ina,...,im) = ay . I] xin; + yi F . (38) 
ij= 
ix] 


Then, INDPGBM*¥ is said to be an interval neutrosophic Dombi power geometric Bonferroni mean 


l —. —— —, — 
Sup (inz, ins) , Sup (inz, ins) is the support degree for nz 
s=1sAz 


from ns, which satisfies the following axioms: (1) Sup (inz,ins) E10; U2) Sup (inz,ins) = Sup (ins,inz); 
(3) Sup (inz,ins) > Sup (ita, inp), if D(inz, ins) < D(ia,iny), in which D(nq,ny) is the distance 


measure between INNs ing and inp defined in Definition 7. 


(INDPGBM) operator. Where T(in:) 


In order to simplify Equation (38), we can describe 


and call A = (Ay,Az,..., 


a, Let Tl) 


(141 (iz)) 


x 


z=] 


Then Equation (38) can be written as follows: 


(39) 


A,)' is the power weight vector (PWV), such that A, > 0, : Az =1. 


z=1 


15 of 32 


Symmetry 2018, 10, 459 
ak, 
2-1 
—— = 1 ! =—IA; <=IA; 
INDPGBM* (ity, ina,...,im) =e | I] Gz + yin; ; ' (40) 
i= 
i#j 
= L u L —u L—L ; 
Theorem 11. Let in; = TR;,TR; | iB; ,ID; | [Fa FL] ), (i = 1,2,...,1) be a group of INNs. 


Then the result obtained from Equation (38) is expressed as 


INDPGBM*» (im iM2 554M ) 


if mi Th Th i 
TRi TR; TRi TR; 
] 1 | 1) 
j | 7 / | 7 
/ / / / 
2 / 1 | ; Pr / U / . 
ys fia| etx vy |y y . + yO J oY fit Balen XY a + y : 
xty f/f |i | De t D. i aty  f \ig | _ i a i" (41) 
/ os / IA, = L IA i =, / i | IA, = IA, =F 
/ / 1- ID; 1-1D; / / 1-1D: 1-1D; 
i] 1 
/ / / 
Pal po ite, of x Pol ae x y 
a / ‘ = ‘ 
fit a* \ y 1 —+ =a yo fit a y | =aoart a 
/ i 1-FL; 1- FL; / 1- FL, 1- FL; 
L u ——T ; 
Theorem 12. (Idempotency) Let inj TR, TR, ,|ID;,1D; |,|FL;,FL;|),@G = 
1,2,...,1) be a group of INNs, if all inj(i i=1,2,...,1) are equal, that is in; = in = 


Qi 
L u L 
ire’ TR al i , 1D |. [Fl ,FL Y= Bae te 


INDPGBM*Y (im, iN, .-+, in) = in. (42) 


Theorem 13. (Commutativity) Assume that in’ , is any permutation of iny (u =1,2,...,1), then 


INDPGBM*! (in, in'y,... in’) — INDPGBM*! (im, ino, ... in) (43) 
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= LE u L =u L =. : 
Theorem 14. (Boundedness) Let inj; = (\TR)LTR | 1D; 10; I [FE:FE;] )( =: D2) kab) 


be a group of INNs, and in = hax TR ,TR; J min 5; , 1D; J min [Fa FE, |) = 
i= i= i= 
mn TR,;,TR; max ID; , 1D; 7max FL;,FL; | ), then 
1= 1= 1= 
= —= = = + 
ino < INDPGBM(iny,i,...,i71) <in . (44) 
: E u L u Ee. ; 
Definition 17. Let in; = ( | TR; ,TR; i 1D; ,1D; | FE; FL; | , (i = 1,2,...,1) be a group of INNs, 
then the INWDPGBM operator is defined as 
_ a P=] 
Iw; (T (inj) +1) leo; (T (inj) +1) 
i= = 1 1 = E we(rie)+1) — E we(tGie)+1) ‘e 
INWDPGBM (iny,ing,...,in) a I] xin; + yin; . (45) 
j= 
if-j 


— L u L =u L—=tL 
Theorem 15. Let in; = (|FR.TR | 1; ,1D; | [FE;,FE;] ). (i = 1,2,...,1) be a group of INNs. 


Then the aggregated result from Equation (45) is expressed as 


INWDPGBM** (in Tagen sini ) 


If 1 
| | / yy / / 
| / | | | / 
(Pat op | teat x ) / | (e- fla of x y 
={ |i fie] y = a -H)) oy fief xy = -+—__? : 
| } es / B | Iw, (rin) =a ws [P01 ea / re / B ffm (rin) (\-7n | ig(Tin +) “f 
\ / / | Yv,(riin)+1)\ ik, Dw, (rin) 4] iR, / Yw.(riin+1)\ mm | Yw,(ri)+1) ik, 
/ | Re: a ot al I 
/ / J / / / 
fof ri en 1 
i ; | | | / ; | | {|| (46) 
P-1 / . | x y } P-l [ | Xx y | 
yo |iey}—x y 2 a Joy fife oy \ = — -HII fh 
| & | 2 [| tu(rinyst) (BY, (rdnysi) ( HY || l = | > w(rinyst) ( ape Ym, (reiny +l) ( HP Y || 
/ | i 7 fr | 7 = j\-w / f 7 = ed i = i =| ||| 
| | » {Tin.)+1)\1- 1D; » {Tin.) +1) 1-1, / » (Pin) +1) 1D: »(T(in.)+1)\1=1D; } \]) 
i / | Le. inz)+ Ym inz) + J) / [ae inz)+ de | in:) + i} J) 
( i 
| i If Wl | f a 
/ / | III | / 
/ } | | | 
f |is(Pats i y i - + , HI i (Eh if y i 7 + u | 
/ Re / a | b(n) +) a i win) +1) iL | | } x+y / ia f (rin) i Zz 7 (rin) a j | 
| | | Wy Tin.) +1 iL} © ng ya\-Foi UE Af | / ¥ Tin.) +1 iE} Sw Tin,) +1 \-FL; J | 
/ | / Pa inz)+ )\ / Yori) ))} | / / | Ye, in.) +1) Xm (in: +1) ) ))) 
\ 


5. MADM Approach Based on the Developed Aggregation Operator 

In this section, based upon the developed INWDPBM and INWDPGBM operators, we will 
propose a novel MADM method, which is defined as follows. 

Assume that ina MADM problem, we need to evaluate u alternatives M= { Mu, Mb, oe mM, 


with respect to v attributes C= {C1,G,. ueay, and the importance degree of the attributes 
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Uv 
is represented by @ = (@1,@2,...,@p)', satisfying the condition @, € [0,1], Yo @,=1. 


The decision matrix for this decision problem is denoted by D = [an] , where den —— 
mxn 


L U L U L u Aa 
(FR FR ; [TD 7D, j [FE p.FEsa| ) is an INN for the alternative My with respect to the 


attribute om (¢g =1,2,...,u;h =1,2,...,v). Then the main purpose is to rank the alternative and 
select the best alternative. 

In the following, we will use the proposed INWDPBM and INWDPGBM operators to solve this 
MADM problem, and the detailed decision steps are shown as follows: 


Step1. Standardize the attribute values. Normally, in real problems, the attributes are of two types, 
(1) cost type, (2) benefit type. To get right result, it is necessary to change cost type of attribute 
values to benefit type using the following formula: 


~ =—L =U u if L u 
dyn = ( |FlgusFlgn|+ |1— Doh 1 — IDga|, |TReue TRen| )- (47) 
Step 2. Calculate the supports 
Supp (dyn dg) =1-D(dgudgt), (g=1,2,.--,uh, 1 =1,2,...,0), (48) 


where, D (den, dq) is the distance measure defined in Equation (10). 


Step 3. Calculate T (den) 


T (dyn) = \_ Supp (den det), (g = 1,2,.-+,jh,1 =1,2,...,0). (49) 
Pi 
Léh 


Step4. Aggregate all the attribute values ds (h = 1,2,...,v) to the comprehensive value Ry, by 
using INWDPBM or INWDPGBM operators shown as follows. 


Rg = INWDPBM (dgi,dgo,..- go) (50) 


or 


Rg = INWDPGBM (det, da, + dgo): (61) 


Step 5. Determine the score values, accuracy values of Re(g = 1,2,...,u), using Definition 5. 

Step 6. Rank all the alternatives according to their score and accuracy values, and select the best 
alternative using Definition 6. 

Step 7. End. 


This decision steps are also described in Figure 1. 
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Read the decision matrix D 


Figure 1. Flow chart for developed approach. 


6. Illustrative Example 


In this part, an example adapted from [42] is used to illustrate the application and effectiveness of 
the developed method in MADM problem. 

An investment company wants to invest a sum of money in the best option. The company 
must invest a sum of money in the following four possible companies (alternatives): (1) car company 
My; (2) food company My; (3) Computer company M3; (4) An arm company My, and the attributes 
under consideration are (1) risk analysis C;; (2) growth analysis C2; (3) environmental impact analysis 
C3. The importance degree of the attributes is © = (0.35, 0.4, 0.25)". The four possible alternatives 
Mg(g = 1,2,3,4) are evaluated with respect to the above attributes C;,( = 1,2,3) by the form of INN, 
and the IN decision matrix D is listed in Table 1. The purpose of this decision-making problem is to 
rank the alternatives. 


Table 1. The IN decision matrix D. 


Alternatives/Attributes Cy C> Cs 
My ({0.4, 0.5], (0.2, 0.3], (0.3, 0.4]) ([0.4, 0.6], [0.1, 0.3], [0.2, 0.4]) ({0.7, 0.9], (0.7, 0.8], (0.4, 0.5]) 
Mp ([0.6, 0.8], [0.1, 0.2], [0.1, 0.2]) ([0.6, 0.7], [0.15, 0.25], [0.2, 0.3]) ([0.3, 0.6], [0.2, 0.3], [0.8, 0.9]) 
Mz ({0.3, 0.6], [0.2, 0.3], (0.3, 0.4]) ([0.5, 0.6], [0.2, 0.3], [0.3, 0.4]) ([0.4, 0.5], [0.2, 0.4], [0.7, 0.9]) 
Mg ({0.7, 0.8], [0.01, 0.1], [0.2, 0.3]) ([0.6, 0.7], [0.1, 0.2], [0.3, 0.4]) ({0.4, 0.6], [0.5, 0.6], [0.8, 0.9]) 


6.1. The Decision-Making Steps 


Step1. Since C1, C are of benefit type, and C3 is of cost type. So, C3 will be changed into benefit 
type using Equation (47). So, the normalize decision matrix D is given in Table 2. 
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Alternatives/Attributes 


Table 2. The Normalize IN decision matrix D. 


C1 


Ce 
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Cs 


My ({0.4, 0.5], [0.2, 0.3], (0.3, 0.4]) (0.4, 0.6], (0.1,0.3], [0.2,0.4])  ([0.4,0.5], [0.2, 0.3], [0.7,0.9]) 
Mp ({0.6, 0.8], [0.1,0.2],[0.1,0.2])  ([0.6, 0.7], (0.15, 0.25], [0.2,0.3]) _({0.8, 0.9], (0.6, 0.7], (0.3, 0.6]) 
Ms ({0.3, 0.6], [0.2, 0.3], (0.3, 0.4]) ({0.5,0.6], (0.2,0.3], [0.3,0.4])  ([0.7, 0.9], [0.6, 0.8], [0.4, 0.5]) 
Mg ([0.7,0.8], [0.01,0.1],[0.2,0.3]) (0.6, 0.7], [0.1,0.2], (0.3,0.4]) [0.80.9], (0.4, 0.5], [0.4, 0.6]) 


Step 2. Determine the supports Supp (Gen, det), (¢ = 1,2,3,4;h,1 = 1,2,3) by Equation (48) (for 


simplicity we denote Supp (dei. det) with S® 


aii el) we have 


Sag = Sid = O50 Si5 = Sis = 0.800; S13 = Sai = 0.85;S}\1 = Bat 7 0,933; 555 43 7 Siam = OFT Sigg = Sri = 0.683; 
SH = Sirs = 0.967; 5,515 = Su = 0.733; S13 = Sia = 0.700; Si\.1> = Sai = 0,902; S15 = Sig = 0.783; S13 = Sits = 0.752; 


Step3. Determine T (den); (¢ = 1,2,3,4;h = 1,2,3) by Equation (49), and we get 


P1800, T= 1.750, 1 = 1600, Fe 161712, = 1600/7 p= A080; 
Tey S067, te = 1.700, 1, = 1A88 1. 160321), — 08d, Ly = 1.905; 


Step4. (a) | Determine the comprehensive value of every alternative using the INWDPBM 


operator, that is, Equation (50) (Assume that x = y = 1;7 = 3), we have 
Ry = ([0.3974, 0.5195], (0.1823, 0.3023], [0.3353, 0.4796]) ; 


Ry = (0.6457, 0.7954], [0.1700, 0.2885], (0.2044, 0.3265]); 
R3 = (0.4846, 0.6503], [0.2556, 0.3711], [0.3376, 0.4394)); 


Rg = (0.6938, 0.7953], (0.1062, 0.2154], [0.3069, 0.4278). 


(b) | Determine the comprehensive value of every alternative using the INWDPGBM 
operator, that is Equation (51), (Assume that x = y = 1; 7 = 3), we have 


x 


Ry, = (0.4026, 0.5381], [0.1570, 0.2977], (0.2998, 0.4520]); 


Ry = (0.6654, 0.8193], (0.1558, 0.2686], [0.1836, 0.3035]); 
R3 = ((0.5159, 0.6732], (0.2366, 0.3473], [0.3265, 0.4279]); 


Rg = ([0.5159, 0.8193], (0.0938, 0.1952], [0.2862, 0.4037]). 

Step 5. (a) Determine the score values of R¢(g = 1,2,3,4) by Definition 5, we have 
S(R1) = 1.8087, S(Rz) = 2.2259, S(R3) = 1.8656, S(R4) = 2.2164; 

(b) Determine the score values of Rg(g = 1,2,3,4) by Definition 5, we have 
S(R1) = 1.8671, S(R2) = 2.2866, S(R3) = 1.9254, S(R4) = 2.1781; 


Step 6. (a) According to their score and accuracy values, by using Definition 6, the ranking order 
is My > My > M3 > My. So the best alternative is M2, while the worst alternative 


is Mj. 
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(b) According to their score and accuracy values, by using Definition 6, the ranking order 
is My > M4 > M3 > My. So the best alternative is M2, while the worst alternative 
is My. 


So, by using INWDPBM or INWDPGBM operators, the best alternative is Mp, while the worst 
alternative is My. 


6.2. Effect of Parameters y, x and y on Ranking Result of this Example 


In order to show the effect of the parameters x and y on the ranking result of this example, we set 
different parameter values for x and y, and = 3 is fixed, to show the ranking results of this example. 
The ranking results are given in Table 3. 

As we know from Tables 3 and 4, the score values and ranking order are different for different 
values of the parameters x and y, when we use INWDPBM operator and INWDPGBM operator. 
We can see from Tables 3 and 4, when the parameter values x = 1 or 0 and y = 0 or 1, the best choice is 
Mg and the worst one is My. In simple words, when the interrelationship among attributes are not 
considered, the best choice is M4 and the worst one is M,. On the other hand, when different values 
for the parameters x and y are utilized, for INWPBM and INWDPGBM operators, the ranking result is 
changed. That is, from Table 4, we can see that when the parameter values x = 1, y = 1, the ranking 
results are changed as the one obtained for x = 1 or 0 and y = 0 or 1. In this case the best alternative is 
M> while the worst alternative remains the same. 


Table 3. Ranking orders of decision result using different values for x and y for INWDPBM. 


Parameter Values INWDPBM Operator Ranking Orders 
T= Ly =O7=3 Sie) = 20606 8(R,) = 24202, Mla > Me > Ns > Mh 
petyesas SRY SISSS SIR)=2288 ay 5 sy > 
cetyetinns SRST SEY 2IE ty sy > 
pesyeinyas ISIS S)=225% ss fa> 
etyeiiras — SEITUB, SR) 226 5 p> > 
T= 1y=47=3 SR) = pa796,s(Rs) = 22552, > Ma > Ma > Mh 
esyereyas — SRITUR SR) 225% ys figs > Bh 


Table 4. Ranking orders of decision result using different values for x and y for INWDPGBM. 


Parameter Values INWDPGBM Operator Ranking Orders 
errs S(R,) = 1.5032, S(Ry) = 1.7934, ar Te ae 
aie eee aad S(Rg) = 1.5136, $(R4) = 1.8037; Bg > Dg tg > Ws 
ee S(Ry) = 1.8220, S(Ry) = 2.2256, oe a 
sae ae Boar S(R3) = 1.8717, S(R4) = 2.1140; td Mei 
se eee S(R1) = 1.8539, S(Ro) = 2.2686, a ee Te 
pa S(R3) = 1.9094, S(R4) = 2.1584; Mee Ma oil os 


ae Ro) = 2.274 ee 
x=5,y=10,7=3 5( Bees) ole 2. My > Mg > M3 > My. 
oe S(R1) = 1.7814, S(Rp) = 2.1710, a ee eer 
se ia S(R3) = 1.8248, S(Ry) = 2.0632; Map a> Mae Mis 
ee S(R1) = 1.8087, S(R) = 2.2259, a ae 
cS N= S(R3) = 1.8656, S(R4) = 2.2164; ae Mae Me: 


WN 
— 


= 2.2866, re ee eee 
Ma > M2 > M3 > Mj. 


ie) 


3) = 1.9254, S(R4 
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From Tables 3 and 4, we can observe that when the values of the parameter increase, the score 
values obtained using INWDPBM decrease. While using the INWDPGBM operator, the score values 
increase but the best choice is Mz for x = y > 1. 

From Table 5, we can see that different ranking orders are obtained for different values of y. When 
y = 0.5 and y = 2, the best choice is M4 by the INWPBM operator; when we use the INWPGBM 
operator, it is M2. Similarly, for other values of y > 2, the best choice is Mz while the worst is Mj. 


Table 5. Ranking orders of decision result using different values for . 


Parameter Values INWDPBM Operator INWDPGBM Operator Ranking Orders 

gaia —0e S(R1) = 1.6662, S(Rz) = 2.1025, S(R1) = 1.7870, S(Rz) = 2.2347, My > Mz > M3 > My. 
i S(R3) = 1.7606, S(Rq) = 2.1972; S(R3) = 1.9103, S(Ry) =2.1812; = My > My > M3 > Mh. 
ae S(R) = 1.7783, S(Rz) = 2.2015, — S(Ry) = 1.8491, S(Rz) = 2.2786, = My > Mp > M3 > Mh. 
Mee S(R3) = 1.8408, S(R4) = 2.2091; S(R3) = 1.9213, S(R4) = 2.1799; = My > My > M3 > Mj. 
cyt S(R1) = 1.8229, S(Rz) = 2.2363, — $(Ry) = 1.8740, S(Rp) = 2.2856, = My > My > M3 > My. 
Co is5 S(R3) = 1.8803, S(Ry) = 2.2219; S(R3) = 1.9275,S(Rq) =2.1751; = My > My > M3 > My. 
pang ae S(R1) = 1.8375, S(Rz) = 2.2455, S(Ry) = 1.8747, S(Rz) = 2.2763, = Mp > My > M3 > My. 
ra S(R3) = 1.9037, S(Ra) = 2.2315; S(R3) = 1.9331, S(R4) = 2.1669; = My > My > M3 > My. 
yy =i is S(Ry) = 1.8418, (Ro) = 2.2477, S(Ry) = 1.8701, S(Ry) = 2.2698, My > Ma > M3 > Mh. 
S(R3) = 1.9160, S(Rq) = 2.2365; S(R3) = 1.9373, S(R4) = 2.1622; My > My > M3 > Mj. 

ee ee ee S(R1) = 1.8447, S(Rz) = 2.2488, S(R1) = 1.8642, S(Rz) = 2.2637, = My > My > M3 > My. 
Bee S(R3) = 1.9270, S(Ry) = 2.2409; S(R3) = 1.9414, S(Rq) = 2.1582; = My > My > M3 > My. 
eye S(R) = 1.8460, S(Rz) = 2.2492, S(Ry) = 1.8608, S(Rz) = 2.2604, = My > My > M3 > My. 
Ci S(R3) = 1.9328, S(Ra) = 2.2432; S(R3) = 1.9435, S(R4) = 2.1562; My >My > M3 > Mj. 


6.3. Comparing with the Other Methods 


To illustrate the advantages and effectiveness of the developed method in this article, we solve 
the above example by four existing MADM methods, including IN weighted averaging operator, IN 
weighted geometric operator [12], the similarity measure defined by Ye [15], Muirhead mean operators 
developed by Liu et al. [42], IN power aggregation operator developed by Liu et al. [37]. 

From Table 6, we can see that the ranking orders are the same as the ones produced by the existing 
aggregation operators when the parameter values x = 1,y = 0,y = 3, but the ranking orders are 
different when the interrelationship among attributes are considered. That is why the developed 
method based on the proposed aggregation operators is more flexible due the parameter and practical 
as it can consider the interrelationship among input arguments. 


Table 6. Ranking order of the alternatives using different aggregation operators. 


Aggregation Operator Parameter Score Values Ranking Order 
INWA operator [12] No . R) _ cet a Se Mg > Mp > Mz > My. 
INWGA operator [12] No : a ~ ae ao amit Mg > My > M3 > My. 

Peer ee * . * zs 

ciminymeme gg} ORE EWRRS = O88 iy 
—— S184, SRB iy iy 
INWMM operator [42] p it 1) : Ry pees Se = en Mg > Mp > M3 > M,. 
INWDMM operator [42] 3 rei 1) : 2 . aoe : ios - Hore My > My > My > My. 
Pope NWOPEM SR) = 10S) 2A ay sy > 
Peper NWO SIDA SUED IME ys iy i> 
INDE porns Yes ARIS) = 2, iy i> 
DuWop ins SRB IR Ry > 
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From the above comparative analysis, we can know the proposed method has the following 
advantages, that is, it can consider the interrelationship among the input arguments and can relieve 
the effect of the awkward data by PWV at the same time, and it can permit more precise ranking order 
than the existing methods. The proposed method can take the advantages of PA operator and BM 
operator concurrently, these factors makes it a little complex in calculations. 

The score values and ranking orders by these methods are shown in Table 6. 


7. Conclusions 


The PBM operator can take the advantage of PA operator, which can eliminate the impact of 
awkward data given by the predisposed DMs, and BM operator, which can consider the correlation 
between two attributes. The Dombi operations of TN and TCN proposed by Dombi have the edge 
of good flexibility with general parameter. In this article, we combined PBM with Dombi operation 
and proposed some aggregation operators to aggregate INNs. Firstly, we defined some operational 
laws for INSs based on Dombi TN and TCN and discussed some properties of these operations. 
Secondly, we extended PBM operator based on Dombi operations to introduce INDPBM operator, 
INWDPBM operator, INDPGBM operator, INWDPGBM operator and discussed some properties of 
these aggregation operators. The developed aggregation operators have the edge that they can take 
the correlation among the attributes by BM operator, and can also remove the effect of awkward data 
by PA operator at the same and due to general parameter, so they are more flexible in the aggregation 
process. Further, we developed a novel MADM method based on developed aggregation operators 
to deal with interval neutrosophic information. Finally, an illustrative example is used to show the 
effectiveness and practicality of the proposed MADM method and comparison were made with the 
existing methods. The proposed aggregation operators are very useful to solve MADM problems. 

In future research, we shall define some distinct aggregation operators for SVHFSs, INHFSs, 
double valued neutrosophic sets and so on based on Dombi operations and apply them to MAGDM 
and MADM. 
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Appendix A. Basic Concept of PBM Operator 
Definition A1. For positive real numbers (,(h = 1,2,...,1) and x,y > 0 the aggregation mapping [54] 


x y 
I 
PBM*" (01, 92,---, 01) = — > eo x AM ei) +1) (Al) 
i=1,j=1 X (T(o) + 1) X (T(0) + 1) 
if-j 


is said to be power Bonferroni mean (PBM) mean operator. 
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Definition A2. For positive real numbers ¢,(h = 1,2,...,1) and x,y > 0 the aggregation mapping [54] 


(T(9)+1) 1(T(9)+1) 
x 1 ! a (T(~0)+1) a (T(~0)+1) 
i mac ere | OO +yo% (A2) 
i=1j=1 
Fj 


is said to be power geometric Bonferroni mean (PGBM) mean operator. 
l 
In Definitions Al and A2,T(9;) = YL supp(, 9;), and supp(9;, @;) is the SPD for ¢; from 9; 
j=l jAi 
satisfying the axioms as; 


(1) sup(~i, pj) = 1 — D( i, 9;), 80 sup (wi, 9;) € [0, 1]; 
(2) sup(@i, 9j) = sup(w;, 8i); 
(3) sup(@i, 9) = sup( Wc. Pa), if |i — %y| < |e — wal: 


where D((j, 9;) is the distance measure among 9; and 7j. 


Appendix B. Proof of Theorem 6 


Proof. Since 


and 
IA,inj =(_|1- : p= : a : ple : =} : pH : 7 
1+] 1A, TR 1+] 1A TR I+] 1A, 1D, I+] 1A, =D 1+] 1A, 1 Fly 1+] 1A, 1 Ly 
1-TR; 1-TR; 1D; 1D; FL; FL; 
Let 
—zZ. == ——f —— 6) =——f = =U =f = =U 
TRi TR: 1-IDi 1-IDi 1-FLi 1-FLi TR; TR; 
a, = T?7i — sc a A 4, = =U ?Si—~ ==1 sh, = ee hs yy = =v”? 
1-TR: 1-TR; ID; ID; FLi FLi 1-TR; 1-TR; 
L =U —f, ——=$ 
«LEWD;  1HID) | IS FL; _1-FLy 
==, ° ° i —— oJ 
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then, we have 


| — nr eee | Poe 282 
1+(JA,)r¢, 1+(IA,)rd, | |1+(UA;)7 8, 1+(IA,)7 4, 


2: 


a ela 1 ji 1 1 1 1 1 


1 * 1 2 1 2 1 ?, 1 


1+(A,)7a, 1+(1A,)’b, L1+(JA,)’c, 1+(JA,)'d, | 1+(IA,)’g, 1+(IA,)rh, 


=i 1 1 1 1 i : 
(1A,in) = a a — | t T Tl T 7 |h 1 as 1 1 : 
14a7/UA)ia 1437/aF, I ee 1ex?/IN)Fd,, | 1+a7/Ua)ie, 1427/Un,)r A 


Coe 1 1 \, 1 4 1 P 1 4 1 


ery (lose 149° /(to,)Fo, [ ey airs 1+y|(to,Fra,, | ixy’/ (tea), 1ey'/ (t0,)Fh 


Moreover, we have 


(iA,in) ®, (IA,in;) = : 7 : z}I- } ple : mae 
1+(x/IA,a! + y/IA,a) ¥ 1+(x/LAbF + y/LA bi) 1+(x/IAje! + y/IAjch)\¥ — 14(x/tA,d? + y/1A,d?) 


1 1 
1- vl +} )s 
1+(a/IA,g/ +y/IAjg7)) 1+(a/IA,h + y/LA,ni } 


and 
rf 


Y(/A,in) ®, (1A ins) 


j=l 
4 


y 


i,j= 
inj 
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1 
7 7\y 
L L / 
=( |1- f |i+ : 7 ffi A : jl- Yd d+ : 7 ffl u ; 
if 1+(x/1A,a7 +y/Lo,a )* / 1+(x/IA,q +y/IA a’ r if + (x/1A,bF +y/IA,b; | 1+(x/1A,by + y/IA bY Yr 


T1 


/ 7 I y 
J I / 1 : I / rl 
Yo J1+] S| t-14+ 7|/}i- ; » Yl 1+] Sy] t-14+ ~| / jie ; 
if 1+(x/IAjeF +y/IA,c)) }/ 1+(x/IA,c? +y/LA,c} )r ifr 1+(x/IA,d? +y/lo,di r / 1+(x/IA,d? +y/IA,d? 
J 
I / r / 1 
Ye i+} ¥} | 1-14 7] / fle ; ,Y {i+} SY] i-i4+ ~|/ l= ; 
iW 1+(x/IA,g7 +y/eAjg))? J) (1+ (x/IA,g7 +y/TAjg3 )P ey + (x/lA,hy +y/IA,hy PJ) (1+ (x/LA hy +y/lohy 


x ; c 
l x y 7 l x y - L x y 
=( |1- 1 14 1 : ,l- 1 14 1 - s| 1 1+ 1 + 
pa lAja;  LAja; x IA,by LA ,b; x lAic?  1A,c; 
i#j j 


x 


ij 


7 7 7 
L x y L x y Ll x y 
VY fi+) 1 Ae a bes pee te ya 
4:7 | 1Ad? "1A d? 47] | 1A.9”  1N.97 SIA IA 
J feat peg i, j= ii joj a et ii 


i#j i#j 


So, we can have 


eats 
Pal 


> (1a,in) ®,(IA,in)). 


i,jal 
izj 


1 
1 ryy \r 
L r l 
fais Meet 4) a 11S da) yi e 
l-l in| \lAjay LA,ay a eAjay eA ay 


1 
1 Ly \r 
L i“ l Y 
a a he ae | | Pe eh if |e ; 
-l ija/ \ LA? LA,bF ial EADS LA be 
i#j i#j 


Y 
1 ! x y 
1 14 1-1/ 14 1 + 
P-l 2, lcs az] 


i#j 


ta 
L x y 
1l/ 14 1 + 
Dy lige IA jd; ) 


i#j 
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A, 
1 1\)" \r 
1 1 a 
1 
te) ae) = ye af | gs 
=i imal \lAjey Ae; imal (lAjey Ae; 
inj , iZj ; 
a 
Es ayy \r 
1 d x y ; u x y J 
1 1+] 1-1 14+) )01 _ 4 Voit 1 : 
P=l 4a/ IAW IA Rh Saf \iAke  1A,h? 
i#j i#j 3 
1 Nf 1 
y ; y 1 2 | ‘ . 
=(|/1- 1 1-1 ; 1+ 1/ 
4] Dt oe IA.a 4] aod lite al Fd Es aA 
itj inj ey 
1 1 1 
Y 1 | . r ie 
124 : Lol : | y |i4|— x aa Ole ; ‘| 
"Pag 1A,d? "Ad Pai iA.g? IA," Pag IA, hi ey 
Then, 
psi 
1 i p F x+y 
pay (Nm) ®, (1A jin,) 
ay 
“ fa on eee ||| Pot lf Safe ge Mae af sal Soha fea 
~ yrs ne / by feel H r ery. / By fameal al i Lye / s) ts mall i (A3) 
i | y f P-l / ii x y i Pai / iif x y fi 
1- If 14; — V 1- If 1+) ——x ]/ 1 7 1- 1 1 V — 
ie y | mlb leear z,]| : ey. / Es ra “| js my / Bs a IA | 
Now, put 
TR TR 1-ID 1-1D, 1 er 12 FE TR; TR 
aj = =, = nn i = 4) = 8 = ar i = fo bj = — tr 
1-TR; 1-TR; : b é FL; 1-TR; 1-TR; 
L — =— 
1-ID; 1-ID; 1-FL; 1—-FL; 
c= = —-, 8} — — 
J] =r 7") =U 7S] an ae a4 | —U 
j i j EL; 
in Equation (A3), we can get 
Bh 
xty 
1 1 ee fahsg 
aa D(A) ®, (1A, in,) 
bia 


iAj 
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Appendix C. Proof of Theorem 7 
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then, there are 


U- U U+ 
FL <FL (n)<FL , 


L- :* L+ U- U U+ L- he L+ U- af U+ L- L L+ 
TR <TR (n,)<TR ,TR <TR (n,)<TR JID <ID(n)<ID ,ID <ID (n)<ID ,FL <FL(n)<FL , 


7 (i) 
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Then there are the following scores 


L U L U L U 
TR +TR ID +ID +1 FL +FL _ 


/ =i | U 
[|e 1-FLi 1-FL; / / - FLi 
/ =r — | / / =v =v 
/ FLi FL; / / / FLi FL; 
| / / 
L- U- L+ U+ Lt U+ 
S TR +TR 4: ID +ID_, 1 FL +FL 
> : 
2 2 2 


Therefore according to the Definition 6, we have 


in < INDPBM(iny, iny,. ..,in1). 


In a similar way, the other part can be proved. That is in < INDPBM (iny,ing,. de im) < 


in . Hence _ 


in < INDPBM(iny, inz,...,im) <in . 
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